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ABSTRACT 

We study the distribution of long Gamma Ray Bursts in the -Epcak — -E-iso and in the ii'peak" 
Fluence planes through an updated sample of 76 bursts, with measured redshift and spectral 
parameters, detected up to September 2007. We confirm the existence of a strong rest frame 
correlation i^poak ct £;0-^4±o.oi (Contrary to previous studies, no sign of evolution with red- 
shift of the i?pcak — ^-iso Correlation (either its slope and normalisation) is found. The 76 
bursts define a strong i^pJ^^j^-Fluence correlation in the observer frame (-Ep^^k °^ -Pboi^ °^) 
with redshifts evenly distributed along this correlation. We study possible instrumental se- 
lection effects in the observer frame i?°J^^]^ -Fluence plane. In particular, we concentrate on 
the minimum peak flux necessary to trigger a given GRB detector (trigger threshold) and 
the minimum fluence a burst must have to determine the value of i^p^ak (spectral analysis 
threshold). We find that the latter dominates in the E'pJ^lj^-Fluence plane over the former Our 
analysis shows, however, that these instrumental selection effects do not dominate for bursts 
detected before the launch of the Swift satellite, while the spectral analysis threshold is the 
dominant truncation effect of the Swift GRB sample (27 out of 76 events). This suggests that 
the i?°J?|]j -Fluence correlation defined by the ^re-Swift sample could be affected by other, still 
not understood, selection effects. Besides we caution about the conclusions on the existence 
of the i?pj?^]^ -Fluence correlation based on our Swift sample alone 
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1 INTRODUCTION 

One of the properties of long Gamma Ray Bursts (GRBs) that re- 
mains mysterious, but potentially fundamental for understanding 
their physics, is the observed correlation between the bolometric 
energy _Eiso emitted during their prompt emission and the the peak 
of the spectrum -Epcak in a vFv plot. In the observer frame a corre- 
lation between the total fluence and iJpcak was found by Lloyd, 
Petrosian & Mallozzi (2000, LPMOO hereafter) with a sample 
of BATSE bursts without measured redshifts. In their pioneering 
work, LPMOO predicted the existence of an intrinsic Speak oc E^^„ 
correlation with a G [0.47, 0.62]. Later, Amati et al. (2002) indeed 
found such a correlation, based on a sample of 12 GRBs observed 
by the BeppoSAX satellite and with spectroscopically measured 
redshifts. The following updates (e.g. Lamb, Donaghy & Graziani 
2005; Amati 2006; Ghirlanda et al. 2007) confirmed this correla- 
tion and showed that the exponent depends somewhat on the fitting 
method and on the sample under consideration, but is of the order 
ofa~0.5. 



E-mail: giancarlo. ghirlanda @ brera. inaf . it 



The current debate about the -Epcak — Eiso correlation con- 
cerns (a) the very existence of the correlation and the presence of 
outliers; (b) the evolution with redshift of the slope and normalisa- 
tion of the correlation and (c) the presence of selection effects on 
this correlation. 

The real existence of the Amati correlation has been ques- 
tioned by Nakar & Piran (2005) and Band & Preece (2005) who 
considered different samples of GRBs, detected by BATSE, of 
known -Epoak but without redshift determination. By considering 
all possible redshifts, they claimed that a large fraction of GRBs 
were in any case outliers to the original (Amati et al. 2002) Am- 
ati relation. They also claimed that the iJpcak — £iso correlation 
was a boundary of a larger dispersion of points in the rest frame 
Speak — Eiao plane. Ghirlanda et al. (2005), using a large sam- 
ple of 442 GRBs with pseudo-redshifts determined by the lag- 
luminosity relation (Norris, Marani & Bonnell 2000; Norris 2002; 
Band & Preece 2005), found that Speak and Siso strongly correlate. 
The correlation has a normalisation and scatter slightly larger than 
in the original Amati et al. (2002) paper. Although these pseudo- 
redshifts are based on the still uncertain lag-luminosity correla- 
tion, they could be used to show that there is no outlier of the 
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Speak — Eiso correlation within the used sample of 442 GRB. Sim- 
ilar conclusions were reached by Bosnjak et al. (2007), using a dif- 
ferent method. Consider that in all cases two GRBs (GRB 980425 
and GRB 031203) were not included in the fits, being clear out- 
liers of the correlation and also anomalous in many ways, even if 
there have been attempts to make them mainstream by considering 
them off-axis events (Ramirez-Ruiz et al. 2005), or sources whose 
prompt flux is dimmed by some scattering material, or spectral evo- 
lution (Ghisellini et al. 2006). 

Li (2007, hereafter L07), considering a sample of 48 GRBs 
(from Amati 2006; 2007) investigated if the correlation evolves 
with redshift, finding that it does, becoming steeper (i.e. larger a) 
at higher redshifts. 

Very recently, Butler et al. (2007, hereafter B07) claimed that 
the correlation exists, but it is probably the result of a selection ef- 
fect, similar to the correlations often found when considering flux 
limited samples and calculating the correlation between the lumi- 
nosities in different bands. They suggest that by multiplying the 
fluence and the observed peak energy by strong function of redshift 
can induce a correlation in the rest frame. 

In this work we study these issues by updating the sample of 
bursts with known redshifts (§2). We study the evolution of the 
-Bpeak — £'iso Correlation with redshift in §3. In §4 we show the 
existence of a strong SpJ^^j. -Fluence correlation. The GRB sample 
considered is heterogeneous in terms of the instruments that de- 
tected the bursts. This requires a deteailed analysis of the possible 
selection effects in order to understand if the Speak — Siso correla- 
tion is an intrinsic correlation or if it is due to any selection effect. 
We study the instrumental selection effects in the _B°j!|ij -Fluence 
observer frame and discuss their impact on the Speak — Eiso cor- 
relation (§4). The relation between the S°J^|]j -Fluence and the 
Speak — Siso relation is briefly discussed in §5, and in §6 we draw 
our conclusions. 

In this work we focus on the sample of bursts with measured 
redshifts and well defined spectral properties in order to study the 
issue of the redshift evolution (L07) and the selection effects on 
the Speak — Siao Correlation (see also B07). We will study the ob- 
server frame Spg^j^-Fluence plane with larger samples of bursts of 
unknown z in a forthcoming paper (Nava et al. in preparation). 

We use f/o = 70 km s"^ Mpc~\ ^m = 0.3 and ^a = 0.7. 
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Figure 2. Con'elation of the spectral parameters (i?°Jj^j^and a) with red- 
shift (top panels) and between the spectral parameters (bottom left) and the 
fluence (bottom right) for the sample of bursts with known redshift reported 
in Tab.n 



2 THE Speak - Siao PLANE 

The information required to put a GRB on the Speak — Siso plane 
are (i) the redshift, (ii) the spectral parameters and (iii) the fluence. 
These are used to compute the bolometric isotropic energy Siso 
(e.g. see Ghirlanda, Ghisellini & Lazzati 2004) and the rest frame 
peak energy Speak of the uF^ spectrum. Note that the fluence and 
peak energy are of the time integrated spectrum, i.e. integrated over 
the total duration of the burst. 

We have collected all bursts with spectroscopically mea- 
sured redshift and with known spectral properties. 35 GRBs 
were detected by instruments on-board different satellites (CGRO, 
BeppoSAX, Hete-Il) before the Swift satellite was launched in Nov. 
2004 (Gehrels et al. 2004) and 41 events were detected by differ- 
ent satellites (Komis-Wind, Swift, RHESSI, Siizaku, Hete-ll ) since 
the end of 2004 in the so called ''Swift era". Among the latter in 
27 cases the spectral parameters (peak energy and fluence) were 
derived from the analysis of the Swiyz-BAT data. For 19 out of 27 
bursts the spectral parameters were taken from the compilation of 
Cabrera et al. 2007 (hereafter C07) and for the other 8 bursts the 



spectral parameters are from the literaturo We refer to this sample 
of 27 bursts as the "Swift burst sample" and define the sample of all 
the other 49 bursts the non-Swift sample. 

In Tab. [T] we report the redshifts, spectral properties and 
isotropic energy of the sample of the 76 GRBs. This is the most 
updated sample of bursts with published redshift. Speak and flu- 
ence that can be put in the Speak — Siso plane up to September 
2007. For many bursts the time integrated spectrum was fitted with 
the Band model (Band et al. 1993) or with a cutoff power-law. 
There are indications (e.g. Kaneko et al. 2006, hereafter K06) that 
if a sufficiently broad energy range is covered by the available spec- 
trum (e.g. up to few MeV), the time integrated spectrum is prefer- 
entially represented by the Band model. In order to uniformly esti- 
mate Siso, following Ghirlanda et al. (2007), for those bursts with 



Note that almost all the 41 bursts detected since 2005 were detected also 
by Swift. However, only in 27 cases the Svvr/t-BAT spectrum could constrain 
the peak energy (C07). In all the other cases only the detection by the other 
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Figure 1. The 76 GRBs with known redshifts and well measured spectral properties (updated to September 2007) in the rest frame plane. Bursts are divided 
in four redshift bins (as labelled). The lines are the best fit obtained with the least squares method (dotted, short-dashed, long-dashed, dot-dashed from small 
to high redshift values, respectively). The insert represents the slope as a function of the 4 redshift bin. The slope of the correlation defined with the entire 
sample of 76 bursts is also shown (cyan symbol in the insert). 



a spectrum fitted by a cutoff power-law we computed the logarith- 
mic average of E^so derived with this model and the value derived 
with the Band function by fixing the high energy spectral index 
to -2.3 (i.e. the typical value reported in e.g. K06). For the GRBs 
taken from COT, Speak was transformed into linear (and its error 
symmetrized), in order to have the same format of several other ta- 
bles already published in the literature. Note however that for the 
points discussed below, our fits do not weight for the errors on the 
variables. All the burst before 2005 taken from Amati (2006) are 
corrected for the different assumption on Ho (Amati 2006 used 
Ho = 65kms~^Mpc'^). 

In Fig. [T] we show the correlation defined by the 76 bursts 
listed in Tab. [T] Statistical analysis gives a Kendall's tau correla- 
tion coefficient r = 0.68 (ISa significance). The correlation can 



be modelled with a power law. Actually, different fitting proce- 
dures have been adopted in the literature: (a) simple least square 
fit which minimizes the difference between the the data points and 
the model (along the ordinate direction) without wighting for the 
errors of the data points or (b) a fit that weights for the errors 
on both variables (see e.g. Press et al. 1986). The least square 
fit gives log Speak = (-22.59 ± 1.39) + log S°-„''^=^°°^ with 
& x^ — 3.49. The fit obtained by weighting for the errors is 
log Speak = (-26.66 ± 0.46) + log S°-55±ooi with a x^ = 545 
for 74 degrees of freedom. In the latter case, the reduced x^ is 
extremely large and this is due to the sample dispersion (see be- 
low) which is much larger than the statistical errors associated with 
the variables. Such a large x^ also leads to underestimate the er- 
rors of the parameters of the fit. It has been proposed (Reichart et 
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Figure 3. Evolution of the slope of the Spoak ~ ^iso correlation (defined 
as i?iso f^ ^ V with redshift. The circles are the results obtained with 



the 48 GRBs of the L07 sample and the open stars are the results obtained 
with the 76 GRBs of our sample. The fit were performed by weighting for 
the erros on both Speak and E^^a for consistency with the method adopted 
by L07. 



al. 2001, 2005, but see Guidorzi et al. 2005, see also D'agostini 
et al. 2006) that the fit should account for an additional parameter 
which is the sample variance. Amati (2006) performed this fit on his 
GRB sample and found a quite large value of the sample variance 
a = 0,15. Instead, we proposed (Ghirlanda et al. 2004, see also 
Liang & Zhang 2005) that the large dispersion of the data points in 
the i?peak — Eiso plane is due to a third observable, i.e. the jet break 
time. Indeed, by correcting the isotropic energy for the collimation 
angle that can be derived from the jet break time of the optical light 
curve, the dispersion of the Speak — Eiso correlation is greatly re- 
duced (Ghirlanda et al. 2007, Nava et al. 2007). In this work we 
will adopt the least square fitting method (Press et al. 1986 p.499) 
without weighting for the errors on the variables because the sam- 
ple dispersion in the Speak — Siso plane defined by the 76 GRBs 
((T=0.22 - consistent with the values found with smaller GRB sam- 
ples, e,g, Ghirlanda et al. 2004, 2005; Amati et al. 2006, C07) is 
clearly larger than the typical statistical errors associated with ei- 
ther Speak and Eiso({uiogEp) = 0-1 and {uiogEi.J = 0.06). 

With the most updated sample of 76 GRBs reported in Tab.[T] 
we can also verify that no correlation exists between the spectral 
parameters or with the redshift (see Fig. |2ll. In particular, there is 
no correlation between a and z nor between a and S°e|i;(second 
and third panel in Fig. [2}. A different, strong, correlation between 
Qpi of the fit with a single power law and the peak energy of the 
fit with a cutoff-power law has been reported recently (Sakamoto 
et al. 2006; Butler et al. 2007) from the analysis of the 5w!//-BAT 
spectra. Note that in this correlation the two spectral parameters Qpi 
and Speak belong to two different models, i.e. a single power law 
and a cutoff-power law model, respectively. However, it has been 
shown in C07 that this correlation is spurious (Fig. 3 and Sec. 3.2 of 
that paper). Although a correlation indeed appears between the sin- 
gle power law spectral index Opi and Spel^ (of the cutoff-power 
law model - see bottom panel of Fig. 3 in C07), it has no physi- 



cal meaning: it is the result of the attempt of the single power law 
model fit to account for the spectral data above the peak which have 
a smaller uFi, flux. Spectra with lower Speak will have a larger 
fraction of data points above the single power law fit and they will 
cause the fit with the single power law to be softer. This accounts 
for the observed positive correlation. It is then dangerous to use this 
correlation for those cases where Spe|]j is not measured, to derive 
it on the basis of api . 

Recently B07 derived the spectral properties of almost all the 
bursts detected by Swift through a Bayesian method. The method 
proposed in B07 to estimate the peak energy, far outside the energy 
range of BAT/Swift (15-150 keV) is based on the assumption of a 
prior distribution for the observed peak energy which is the (gaus- 
sian) peak energy distribution of BATSE bright bursts (Kaneko et 
al. 2005), for Spelk>300 keV, and a uniform distribution for en- 
ergies below this value. However, as also noted by B07, from the 
comparison of these estimates of Speak with those derived from 
the spectral fits of the Konus-Wind and Suzaku, for the common 
bursts, there could be a bias at high values of E°^ly. estimated from 
Swift data through the Bayesian method and this could influence 
the findings of B07. 

In this paper, we adopt a conservative approach and we con- 
sider only the Swift GRBs for which the estimates of the peak en- 
ergy is based on a standard (i.e. frequentist) fit of the BAT spectra 
of bursts with known redshifts (from C07). 



NO REDSHIFT EVOLUTION OF THE Speak 
CORRELATION 
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The 76 bursts of Tab.[T]are distributed in a large redshift range, up 
to z = 6.3 (for GRB 050904 - Tagliaferri et al. 2005). It is worth 
to investigate if any evolution with z can affect the Speak — Siso 
correlation defined by this sample. Note that the simple evolution of 
the burst energetics (Siso) with z does not necessarily implies the 
evolution of the Speak — Siso correlation. Indeed, if Speak evolves 
in the same way as Siso (i.e. if the link between Siso and Speak is 
physically robust) then we should see no evolution of the slope and 
normalisation of the Speak — Siso correlation with z. 

The possible redshift evolution of the Speak — Siso correlation 
has been investigated by Li 2007 (L07, hereafter). He used the sam- 
ple of 48 GRBs (from Amati 2006, 2007), restricting to long events 
and excluding the peculiar GRB060614 (Gherels et al. 2006) and 
dividing the sample into 4 redshift bins. Comparing the slopes of 
the correlation corresponding to each redshift bin he found that the 
Speak — Siso Correlation becomes steeper for increasing redshift. 
Note that the slope defined by L07 refers to the Siso oc S^^j^ cor- 
relation, i.e. the reverse of that defined in this paper. This is opposite 
to what one naively expects, namely that for larger z we select, on 
average, more energetic bursts with the same Speak, resulting in 
a. flattening of the Speak — Siso slope. Note also from Fig. [T] that 
bursts with different redshifts are distributed along the Speak — Siso 
correlation with no evident segregation along the correlation itself 
(except for the excess of low redshift bursts in the low end of the 
Speak - Siso plane). 

In our study we examine first the possible evolution with red- 
shift of the Speak — Siso Correlation defined as Speak tx Si^o- Our 
sample of Tab.[T]extends the original sample of 48 GRBs, used by 
L07, both in number and redshift. To verify his results, we have 
divided the 76 GRB sample into four redshift bins, chosen to have 
an equal number of 19 bursts per bin. These z bins are very sim- 
ilar to those adopted by L07. We fit the Speak — Siso correlation 
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in each sub-sample through the least square method (see Sec. 2). 
The results are shown in the insert of Fig. [T] The analysis of the 
S'peak — Eisa Correlation in each redshift bin excludes the evolu- 
tion of the slope of the Speak — Eiso correlation with z, given the 
present sample of 76 events. Indeed, the slopes of the iJpcak — £iso 
correlations defined by the four redshift sub-samples are consistent 
among themselves and also with the slope of the correlation defined 
by the entire sample (insert of Fig.[Tll. We have also verified that the 
normalisation of the correlations defined by the four redshift bins 
are consistent. 

We then investigated the possible reasons for the discrepancy 
of our results with respect to L07. We reconstructed the sample 
of 48 GRBs used by LOT (with data taken from Amati et al. 2006, 
2007). Using this sub-sample we find the same results of Li (2007), 
concerning the evolution of the slope and normalization of the 
-Epeak — Eiso Correlation with redshift (circles in Fig. [3](. How- 
ever, when extending this sample to the 76 GRBs of our present 
sample, we do not find any evolution with redshift (star symbols 
in Fig. [3}. Note that for this comparison we used the same fitting 
method of L07, which weights for the errors on both variables, and 
fitted the _Eiso oc -E ^ak correlation. We then conclude that the Li 
(2007) results were affected by too low statistics, and that there is 
no evidence that the iSpcak — Eiso correlation evolves. 



4 THE E: 
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-FLUENCE CORRELATION 



The rest frame Ep^ak — Eiso correlation is defined by multiplying 
the observed peak energy by (1 -I- z) and the observed (bolomet- 
ric) fluence by 'iiTd'l/{l + z), where cJl is the luminosity distance. 
Therefore, in principle, the Speak — Eiso correlation could be only 
apparent, being the result of multiplying Spelk ^^^ the fluence F 
by strong functions containing the redshifts (as argued by B07). 

To study if this is the case, it is necessary to study what region 
of the observational plane Spelj^ -Fluence is accessible to observa- 
tions. Any detector, in fact, can only observe GRBs above a limiting 
fluence and within a limited energy range. 

Other selection effects could be related to the requirement of 
measuring the redshift: for instance (in the absence of X-ray lines) 
we require the burst to be located precisely enough to make possible 
the optical follow up and we also require its afterglow to be visible 
in the optical. If only bright bursts had an optical afterglow, this 
would affect the resulting correlation in the Speak — Eiso plane. 

We focus here on the selection effects introduced by requiring 
that the burst is not only detected, but it must be bright enough to 
allow the determination of the spectrum of its prompt emission. 

Consider first how our 76 GRBs with known redshift are lo- 
cated in the SpeJj^ -Fluence plane, as illustrated in Fig.|4] Different 
symbols correspond to GRBs in different redshift bins. Note that: 

• These bursts define a strong correlation in this plane. The 
Kendall's correlation coefficient is r = 0.53 (9cr significance). The 
least square fit has a slope 0.39±0.04. However, it is very likely 
that this correlation is affected by truncation effects which should 
be considered when computing the correlation strength and when 
recovering its "true" slope (LPMOO). 

• There is no segregation in redshift. We do not find that nearby 
bursts are brighter and bluer, and distant bursts are fainter and red- 
der, as we would expect if there were no intrinsic correlation be- 
tween Speak and Siso. 

The existence of a S°eak -Fluence correlation is not new: it 
was discovered with a relatively small sample of BATSE bursts by 
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Figure 4. i?°^Jj^-boIoraetric fluence plane for the sample of 76 bursts of 
Tab. 1 . We show with different symbols the redshift bins corresponding to 
0.03-0.843 (circles), 0.846-1.48 (stars), 1.49-3.1 (squares) and 3.2-6.29 
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(triangles). The least square fit to the data is E' 

the fluence in the energy range 1 keV - 10 MeV. The scatter of these bursts 

around the best fit has a cr = 0.21. 



LPMOO and subsequently confirmed with the discovery of the class 
of X-Ray Flashes by BeppoSAX and Hete-II (e.g. Lamb, Donaghy 
& Graziani 2005). Sakamoto et al. (2005) showed that a correlation 
between these two observables exists in the Hete-II burst sample 
and extends from the softest S°j!|,^ of a few keV to the few hun- 
dred keV range (e.g. Preece et al. 2000 and K06 for the BATSE 
sample). We will discuss the sample of bursts with no redshift, but 
with measured S°eak' in a forthcoming paper (Nava et al. 2008, in 
preparation). 



4.1 Instrumental biases in the SpeJ^-Fluence plane 

Consider the two regions in the Spelit-Fluence plane of Fig.|4]be- 
low and above the correlation where we have a paucity of points. 
Region 1 is characterized by large fluences and low/moderate 
Speli^; region 2 comprises small fluences, large/moderate Spe^j^. 
In region 1 there are no (or very few) bursts, although there is no 
instrumental effect against their detection. This implies that they 
really do not exist (or, rather, they are very few). In region 2, in- 
stead, the paucity of points in the SpeJj^ -Fluence plane could be 
affected by some instrumental selection effect. 

The detectors (past and present) dedicated to study GRBs in- 
troduce at least two selection effects on the population of bursts 
that they observe: (a) the trigger sensitivity which is the ability to 
detect a burst as a transient event significantly above the average 
background (see e.g. Band 2003; 2006); (b) the "spectral analysis 
threshold" which is the minimum fluence required to constrain the 
peak energy from the spectral analysis. 

For any Speak. ^ minimum photon flux is required to trigger 
the burst. This minimum flux depends on the spectral properties of 
the burst (especially its peak energy - e.g. Band 2003), on the detec- 
tor design and on the total background rate seen by the instrument. 
The ability to detect a GRB for a non-imaging instruments (like 
BATSE, Komis, Hetell-Fregate, SeppoS AX-GRBM), is related to 
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-1.05 [0.21] 


-2.2 [0.1] 


7.5e-5 [0.8e-5] 


40-700 


348 [28] 


1.82e53 [0.23e53] 


1 


990712 


0.433 


-1.88 [0.07] 


-2.48 [0.56 ] 


6.5e-6 [0.3e-6] 


40-700 


93 [15] 


6.72e51 [1.29e51] 


1 


991208 


0.706 






1.6e^[5.0e-6] 


20-10000 


313 [31] 


2.23e53 [0.18e53] 


5 


991216 


1.02 


-1.23 [0.13] 


-2.18 [0.39] 


1.9e^[0.2e^] 


20-2000 


642 [129] 


6.75e53 [0.81e53] 


2 


000131 


4.50 


-0.69 [0.08] 


-2.07 [0.37 ] 


4.2e-5 [0.4e-5] 


20-2000 


714 [142] 


1.84e54 [0.22e54] 


2 


000210 


0.846 






7.6e-5 [5.0e-6] 


20-10000 


753 [26] 


1.49e53 [0.16e53] 


5 


000418 


1.12 






2.6e-6 [4.0e-5] 


20-10000 


284 [21] 


0.91e53 [0.17e53] 


5 


000911 


1.06 


-1.11 [0.12] 


-2.32 [0.41 ] 


2.2e^ [0.2e^] 


15-8000 


1856 [371] 


6.7e53 [1.4e53] 


5 


000926 


2.07 






2.6e-5 [4.e-6] 


20-2000 


310. [20] 


2.7e53 [0.58e53] 


5 


010222 


1.48 






1.4e^ [8.e-6] 


20-2000 


766 [30] 


8.1e53 [0.86e52] 


5 


010921 


0.45 


-1.60 [0.1] 




1.84e-5 [O.le-5] 


2^00 


129. [26] 


I.le52[0.11e52] 


5 


011211'' 


2.140 


-0.84 [0.09] 




2.6e-6 [0.3e-6] 


40-700 


185 [25] 


6.64e52 [1.32e52] 


(7)2 


020124 


3.198 


-0.87 [0.17] 


-2.6 [0.65 ] 


8.1e-6[0.9e-6] 


2^00 


390 [113] 


2.15e53 [0.73e53] 


8 


020405 


0.695 


-0.00 [0.25] 


-1.87 [0.23] 


7.4e-5 [0.7e-5] 


15-2000 


617 [171] 


1.25e53 [0.13e53] 


2 


020813'' 


1.255 


-1.05 [0.11] 




1.0e^[0.1e^] 


30-400 


478 [95] 


6.77e53 [1.00e53] 


3 


0208 19B 


0.41 


-0.90 [0.15] 


-2.0 [0.35 ] 


8.8e-6 [9.e-7] 


2^00 


70. [21] 


6.8e51 [0.17e51] 


(7)5 


020903'' 


0.25 


-1.00 [0.0] 




5.9e-8 [1.4e-8] 


2-10 


3.37 [1.79] 


1.8d49[0.31d49] 


9 


021004* 


2.335 


-1.00 [0.2] 




2.6e-6 [0.6e-6] 


2^00 


267 [117] 


4.09e52 [0.71e52] 


7 


021211 


1.01 


-0.85 [0.09] 


-2.37 [0.42 ] 


2.2e-6 [0.2e-6] 


30-400 


94 [19] 


I.le52[0.13e52] 


2 


030226'' 


1.986 


-0.90 [0.2 ] 




5.6e-6 [0.6e-6] 


2^00 


290 [63] 


6.7e52 [1.20e52] 


7 


030328 


1.520 


-1.14 [0.03] 


-2.1 [0.3 ] 


3.7e-5 [0.14e-5] 


2^00 


328 [35] 


3.61e53 [0.40e53] 


3 


030329 


0.169 


-1.32 [0.02] 


-2.44 [0.08 ] 


1.2e^[0.12e-4] 


3^00 


79 [3] 


1.66e52 [0.20e52] 


3 


030429'' 


2.656 


-1.10 [0.3] 




8.5e-7 [1.4e-7] 


2^00 


128 [37] 


1.73e52 [0.31e52] 


7 


040924'= 


0.859 






2.6e-6 [0.0] 


30-400 


102 [35] 


0.95e52[0.1e52] 


5 


041006'' 


0.716 


-1.37 [0.14] 




2.0e-5 [0.2e-5] 


25-100 


108 [22] 


8.30e52 [1.3e52] 


3 


050126'' 


1.29 


-0.75 [0.44] 




8.55e-7 [1.82e-7] 


15-150 


263 [110] 


7.36e51 [1.60e51] 


10 


050223'*''' 


0.5915 


-1.50 [0.42] 




6.14e-7[0.83e-7] 


15-150 


110 [54] 


1.21e51 [1.77e50] 


10 


050318'' 


1.44 


-1.34 [0.32] 




2.1e-6[0.2e-7] 


15-350 


115 [27] 


2.00e52 [0.31e52] 


11 


050401= 


2.9 


-1.00 [0.0] 


-2.45 ... 


1.93e-5 [0.04e-5] 


20-2000 


501 [117] 


4.1e53 [0.8e53] 


13 


0504 16A 


0.653 


-1.01 [0.0] 


-3.4... 


3.5e-7 [0.3e-7] 


15-150 


28.6 [8.3] 


8.3e50 [2.9e50] 


14 


050505'' 


4.27 


-0.95 [0.31] 




2.58e-6 [3.06e-7] 


15-150 


661 [245] 


1.76e53 [2.61e52] 


10 


050525A'' 


0.606 


-0.01 [0.11] 




2.01e-5 [0.05e-5] 


15-350 


127 [5.5] 


2.89e52 [0.57e52] 


25 


050603 


2.821 


-0.79 [0.06] 


-2.15 [0.09] 


3.41e-5 [0.06e-5] 


20-3000 


1333 [107] 


5.98e53 [0.4e53] 


15 


050803'' 


0.422 


-0.99 [0.37] 




2.08e-6 [2.57e-7] 


15-150 


138 [48] 


1.86e51 [3.99e50] 


10 


050814'' 


5.3 


-0.58 [0.56] 




1.46e-6 [1.16e-7] 


15-150 


339 [47] 


1.12e53 [2.43e52] 


10 


050820A'' 


2.612 


-1.12 [0.14] 




5.27e-5 [1.2e-5] 


20-1000 


1325 [277] 


9.75e53 [0.77e53] 


16 


050904 


6.29 


-1.11 [0.06] 


-2.2 [0.4] 


5.4e-6 [l.e-7] 


15-150 


3178 [1094] 


1.24e54[0.1e54] 


26 


050908'' 


3.344 


-1.26 [0.48] 




4.36e-7 [0.46e-7] 


15-150 


195 [36] 


1.97e52[3.21e51] 


10 


050922C'' 


2.198 


-0.83 [0.26] 




2.6e-6 [0.26e-6] 


30-400 


417 [118] 


4.53e52 [0.78e52] 


17 


051022'' 


0.80 


-1.17 [0.038] . 




2.61e-^ [0.8e-^] 


20-2000 


918 [63] 


5.3e53 [0.5e53] 


27 


051109A'' 


2.346 


-1.25 [0.5] 




4.0e-6 [l.Oe-6] 


20-500 


539 [381] 


7.52e52 [0.88e52] 


18 


060115'' 


3.53 


-1.13 [0.32] 




1.6e-6 [1.07e-7] 


15-150 


288 [47] 


7.38e52 [1.23e52] 


10 


060124'' 


2.297 


-1.48 [0.02] 




2.5e-5 [0.35e-5] 


20-2000 


636 [162] 


4.3e53 [0.34e53] 


19 


060206'' 


4.048 


-1.06 [0.34] 




8.4e-7 [0.4e-7] 


15-150 


381 [98] 


4.68e52 [0.71e52] 


20 


060210'' 


3.91 


-1.12 [0.26] 




6.92e-6 [3.74e-7] 


15-150 


575 [186] 


4.15e53 [5.7e52] 


10 


060218'' 


0.0331 


-1.62 [0.16] 




3.7e-6 [0.37e-6] 


15-150 


4.9 [0.3] 


7.7e49 [1.42e49] 


29(31) 


060223A'' 


4.41 


-1.16 [0.35] 




6.54e-7 [0.52e-7] 


15-150 


339 [63] 


4.29e52[6.64e51] 


10 


060418'' 


1.489 


-1.50 [0.15] 




1.6e-5 [0.16e-5] 


20-1100 


572 [114] 


1.28e53 [0.10e53] 


21 


06051 OB'' 


4.9 


-1.53 [0.19] 




3.86e-6 [2.88e-7] 


15-150 


575 [227] 


3.67e53 [2.87e52] 


10 


060522'' 


5.11 


-0.70 [0.44] 




1.05e-6 [1.04e-7] 


15-150 


427 [79] 


7.77e52 [1.52e52] 


10 
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GRB 


z 


a 


/3 


Fluence 
erg cm~^ 


Range 
keV 


-S-pcak 

keV 


erg 


Ref 


060526-' 


3.21 


-1.1 [0.4] 


-2.2 [0.4] 


4.9e-7 [0.6e-7] 


15-150 


105.25 [21.1] 


2.58e52 [0.26e52] 


26 


060605'' 


3.78 


-1.0 [0.44] 




5.33e-7[1.12e-7] 


15-150 


490 [251] 


2.83e52[4.5e51] 


10 


060607A'' 


3.082 


-1.09 [0.19] 




2.49e-6 [1.58e-7] 


15-150 


575 [200] 


1.09e53 [1.55e52] 


10 


060614 


0.125 


-10 [0] 




2.2e-5 [0.22e-5] 


15-150 


55 [45] 


2.51e51 [1.0e51] 


28 


060707'' 


3.43 


-0.73 [0.4 ] 




1.63e-6[1.13e-7] 


15-150 


302 [42] 


6.62e52[1.39e52] 


10 


060714'' 


2.711 


-1.77 [0.24] 




3.06e-6 [3.96e-7] 


15-150 


234 [109] 


1.34e53 [9.12e51] 


10 


060814 


0.84 


-1.43 [0.16] 




2.69e-5 [0.26e-5] 


20-1000 


473 [155] 


7.01e52 [0.7e52] 


32 


060904B'' 


0.703 


-1.07 [0.37] 




1.48e-6[1.91e-7] 


15-150 


135 [41] 


3.64e51 [7.43e50] 


10 


060906'' 


3.686 


-1.6 [0.31] 




2.38e-6 [1.92e-7] 


15-150 


209 [43] 


1.49e53 [1.56e52] 


10 


060908'' 


2.43 


-0.9 [0.17] 




2.66e-6[l.lle-7] 


15-150 


479 [110] 


7.79e52 [1.35e52] 


10 


060927'' 


5.6 


-0.93 [0.38] 




l.le-6[0.1e-6] 


15-150 


473 [116] 


9.55e52 [1.48e52] 


22 


061007 


1.261 


-0.7 [0.04] 


-2.61 [0.21] 


2.5e^[0.15e-^] 


20-10000 


902 [43] 


8.82e53 [0.98e53] 


23 


061121'' 


1.314 


-1.32 [0.05] 




5.7e-5 [0.4e-5] 


20-5000 


1289 [153] 


2.61e53 [0.3e53] 


24 


061126'' 


1.1588 


-1.06 [0.07] 




3.0e-5 [0.4e-5] 


15-2000 


1337 [410] 


3.0e53 [0.3e53] 


12 


061222B'' 


3.355 


-1.3 [0.37] 




2.24e-6 [1.23e-7] 


15-150 


200 [28] 


1.03e53 [1.6e52] 


10 


070125 


1.547 


-1.1 [0.1] 


-2.08 [0.13] 


1.74e-^[0.17e^] 


20-10000 


934 [148] 


9.3e53 [0.93e53] 


30 



Table 1. " GRB 970508 is included following Amati (2006), but it was excluded in Ghirlanda et al. (2004) and Ghirlanda et al. (2007) because of the 
inconsistency of the BATSE and BeppoSAX spectram (i.e. it is a candidate outHer if the BATSE spectrum is assumed). In the sample of bursts before 2005 
there are some -Epcak and E^^^ values which are different from the recent compilation of Amati (2006). This is because he takes the average when two or 
more instruments (e.g. Konus and Hete) made the spectrum. This is correct when the spectral results are consistent, but we do not adopt the same values if 
the spectral results are inconsistent. * E\^a has been computed as the logarithmic average of the values found with the Band (with (3 fixed at -2.3) and cutoff 
power-law model (see also Ghirlanda et al. 2007, Firmani et al. 2005). '^ We adopt the peak energy and isotropic energy values reported in Amati et al. (2006). 
There are only two GCNs for this burst giving the fluence and peak energy derived by Hete-ll and Konus, which are consistent. We then take the average 
of the two values, a values are not given in the GCN. '' The redshift is known but although it is present in B07 they do not report the redshift. "^ We adopt 
an average of the spectral properties of the two peaks reported by the GCNs. Note that the two spectra are consistent, f This is the burst discussed on our 
webpage (http://www.brera.inaf.it/utenti/gabriele/060526.html I. This burst is an outlier of the Amati relation if the peak energy is much larger. References: 
before 2005: (1) Amati et al. (2002); (2) Ghirlanda et al. (2004, and references therein); (3) Nava et al. (2006 and references therein); (4) Ulanov et al. (2005); 
(5) Amati (2006 and references therein); (6) Firmani et al. (2006); (7) Sakamoto et al. (2005); (8) Atteia et al. (2005); (9) Sakamoto et al. (2004); Since 2005: 
(10) Cabrera et al. (2007). For the error the geometric mean is computed from the extremes of the logarithmic values); (11) Peni et al. (2005); (12) Perley et 
al. (2007, spectral analysis of the BAT+RffSSS/ spectrum); (13) Golenetskii et al. (2005); (14) Sato et al. (2005); (15) Golenetskii et al. (2005a); (16) Cenko 
et al. (2006); (17) Crew et al. (2005); (18) Golenetstkii et al., (2005b); (19) Romano et al. (2006); (20) Palmer et al. (2006); (21) Golenetskii et al. (2006); (22) 
Stamatikos et al. (2006); (23) Golenetskii et al. (2006a); (24) Golenetskii et al. (2006b, Epcak and i?iso are the average of the values found by the RHESSI 
and Konus spectra, while the spectral parameters are of the Konus spectrum); (25) Blustin et al. (2005); (26) Schaefer (2007); (27) Golenetskii et al. (2006c); 
(28) Mangano et al. (2007); (29) Campana et al. (2006); (30) Golenetskii et al. (2007). The burst was also seen by RHESSI (Bellm et al. 2007) whose spectrum 
is inconsistent with that of Konus. However, the two i?poak and the fluences values given in the GCN are inconsistent with the average Epcak reported in the 
same GCN (for RHESSI). For this reason we report and use the Konus spectral data; (31) Ghisellini et al. (2006); (32) Golenetskii et al. (2006d). 



the significance of the burst intensity (i.e. its peak flux) with respect 
to an average background (see e.g. Band et al. 2003). For imagers 
this threshold is modified by the coding aperture of the mask and it 
can be even lower (Band 2006). The trigger sensitivity curves as a 
function of the burst peak spectral energy for different instruments 
has been computed by Band (2003) (see also Band 2006). We show 
them for BATSE, Swift, Hete-II, BeppoSAX in Fig.|5] These curves 
correspond to the sensitivity limit of each instrument to burst with 
a varying _E°j!|j;and are computed by Band (2003) in the peak flux 
- peak energy space. In order to report them on our iSpcai^ -Fluence 
plane we have (i) to assume a typical burst duration to cumpute the 
fluence and (ii) to assume a burst spectrum to convert from the pho- 
ton units to the energy units. To make this conversion we assumed 
a burst with a 1 s duration which represents a lower limit for the de- 
tection sensivity of long GRBs and we also assumed a burst with a 
typical Band spectrum (with a = —1.0 and /3 — —2.3) and with a 
variable Spclk- As discussed by Band 2003, however, these curves 
are only an estimate of the burst detection sensitivity for any given 
instrument because the trigger strongly depends on the burst time 
profile and on the background. In particular we tried to account for 
the unknown burst duration when converting the Band (2003) lim- 
iting curves in the fluence plane of Fig \5\ In the simplest case of 
a burst with a single pulse with a triangular shape, the conversion 



is -F(£pcak) = P(£;pcak)r/2 where i^(£;pcak) and P(£pcak) are 
the fluence and peak flux as a function of the peak energy of the 
burst and T is the burst duration. Under this simple hypothesis, 
our curves for the trigger threshold would move in the plot of Fig. 
[5] proportionally to T/2. However, GRBs light curves are all but 
simple single peaked triangles. A more realistic conversion factor 
between the peak flux and the fluence, which also accounts for the 
typical variation of the flux with time, can be obtained from the ra- 
tio of the fluence and the peak flux F/P for the complete GRBs of 
the BATSE sample (taking into account only the long burst popula- 
tion). We found that the distribution of the ratio F/P peaks at ~6 
and has a tail towards lower values. In Fig.[5]we report the two lim- 
iting curves for the trigger defined with a burst of duration 1 s and 
with the above ratio. Note that the WFC trigger thresholds are on 
the left of the plot and are clearly not affecting the distribution of 
bursts with known redshifts. Moreover, these curves are here plot- 
ted in the very narrow energy range corresponding to the WFC trig- 
ger but their behaviour at higher i5°cak values has been shown to 
be very weak (see Fig. 3 of Band 2003). The WFC thresholds were 
taken from Band (2003) and they are consistent with what reported 
by Vetere et al. (2007JJ. This is contrary to what one could think, 

2 Vetere et al. (2007) report a WFC detection threshold of ~ 4 X 10~^ 
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Figure 5. Observer frame plane of the peak energy and bolometric fluence (1-10* keV). Bursts whose peak energy has been derived by the analysis of the 
BAT data are shown with filled squares. The open symbols are all the other bursts with redshift. The shades regions represent the minimum fluence that a burst 
can have (for different duration of the burst -5s and 20 s for the left and right boundary of each shaded region, respectively) in order to perform the spectral 
analysis and derive constraints on £^°J^^j. ■ For Swift we also report (narrow grey stripe) the limiting curve derived by slightly different criteria (see text for 
details). The single lines represent the trigger sensitivity for Swift (dashed), BATSE (solid), BeppoSAX (dotted), Hete-II (triple dot-dashed). The best fit of 
the sample of Swift bursts (long-dashed line) and non-Swift bursts (dashed line) and of the total sample (solid line) are also reported. These fits, however, do 
not take into account the truncation effects (see text for details). The open red square represents GRB 060218 whose peak energy has been determined by the 
analysis of the SuT/f-XRT data (Campana et al. 2006). 



i.e. that due to the requirement to accurately locate the burst in or- 
der to measure its redshift, it is the WFC trigger threshold to bias 
the GRBs (of the pre-Sw/// sample) with measured redshift in the 
Spclij-Fluence plane (but see also the discussion at the end of this 
section). These results, taken at face values, imply that the distribu- 
tion of GRBs of our heterogeneous sample in the i^p^lj^ -Fluence 



erg cm~^ s~^ in the 2-10 keV range which is in good agreement with the 
curves of Band (2003) when converted in bolometric flux (i.e. in the 1-10'^ 
keV band). 



plane are not affected by the trigger sensitivity of the different istru- 
ments that detected these bursts. 

However, in order to perform the spectral analysis and to con- 
strain any model parameters we also require a minimum number 
of total (time integrated) photons. To derive the threshold fluence 
^min required to perform the spectral analysis, we have simulated 
several spectra with different E'^l^ and fluence values. For this 
simulation we need the detector response and a typical background 
spectrum. These two files are available for BATSE and Swift at 
the respective web sites. Unfortunately, no public response file was 
found for the BeppoSAX-GRBM instrument. For the Swift simu- 
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lations we used the batphasimerr tool of the latest HASOFTv6.3.2 
release (C. Markwardt, private communication). 

We made two sets of simulations, assuming a burst durations 
of 5 and 20 seconds, respectively, with a corresponding different 
amount of background. 

For each simulation we assumed an intrinsic spectrum de- 
scribed by the Band model with typical spectral slopes (see e.g. 
K06) a = —1 and /3 = —2.3, a given peak energy and a given 
bolometric fluence. The resulting simulated spectrum is then fitted 
with a Band model (for BATSE bursts) or with a cutoff-power law 
model (CPL, for Swift bursts) and with a simple power law (PL). 
The analysis returns the values of the spectral parameters and the 
corresponding value of the reduced x^, for both the Band or CPL 
model and for the PL one. We then repeat this procedure 500 times, 
constructing the distributions of the "output" values of -Bpelk (f™' 
the Band and CPL models) and the associated error. We then fit 
the two distributions with a Gaussian to find the average value of 
^peak' the average value of the error with the associated cte value. 
The average value of the "output" -Ep^ak 's consistent with the as- 
sumed -Bpcak value. We then decide if the fitting is returning an 
acceptable value of E°]^ly. adopting two criteria: 



• The error on -Epelkis less than 100% in the majority (97.7%) 
of the 500 simulations. 

• The fit with the Band or CPL models are significantly better 
than the fit with a single power law. This is measured by the id- 
lest, and we set the threshold to 2a (corresponding to a probability 
of 95.45%). 

Note that for Swift bursts we do not ask that an intrinsic spec- 
trum following the Band function can be better fitted with the same 
model, requiring only that we can reconstruct -Ep^ak' even if only 
with a CPL model. If both these conditions are fulfilled, we repeat 
the set of simulations with a smaller fluence and the same -EpJ?ak> 
until one or both the conditions are not met. This defines the min- 
imum fluence (for a given E°^ljJ. This is illustrated in Fig. [5] by 
the shaded regions (they are not lines because they correspond to 
two different duration times). The values reported in Fig. [5] are the 
"output" values of -Bpclk ^rid bolometric fluence, since we want to 
compare them with the GRBs actually detected. Our simulations 
have been performed assuming a peak energy between 30 and 140 
keV for Swift and 50 and 1000 keV for BATSE. For values of Epeak 
outside this range it is difficult to fit a cutoff-power law or Band 
spectrum due to the few spectral bins between the peak energy and 
the end of the energy range. 

In another set of simulations, performed only for Swift bursts, 
we relaxed the above criteria, replacing them by the condition that 
the "output value" of -Bpelk determined by a cut off power law 
model, and its associated Icr error, is contained in the 15-200 keV 
band. In other words, E°]^lj^+aE has to be less then 200 keV, and 
E^pcak^'^B must be larger than 15 keV. For these simulations we 
performed the same spectral fitting procedure described in Section 
2 of C07 (in particular a logarithmic spacing of -Bpeak was used). 
The resulting limiting curve is shown by the grey narrow stripe in 
Fig. |5] (corresponding to burst durations of 5 seconds). 

Note that the simulations assume a burst observed on-axis for 
both Swift and BATSE. While in the case of BATSE the LAD sen- 
sitivity does not strongly change for off-axis incidence, in the case 
of Swift-BAT the sensitivity strongly depends on the photon inci- 
dence angle. For this reason we also considered for Swift the case 



of a 60 deg off-axis burst (light shaded region in Fig. 151 I We also 
tested the dependence of the simulations from the assumed shape 
of the spectrum: we repeated the simulations with a = —0.5 (i.e. 
in the tail of the distribution of this parameter for the BATSE bright 
burst population - K06) and found that the limiting curves in Fig. 
[5]move to the left by a factor 2 on average. 

To understand the shape of the limiting region in Fig.[5]we can 
make a simple argument. It is reasonable to assume that the deter- 
mination of -Epoak (within a given confidence level) requires a min- 
imum number of photons A'^min around the peak of the spectrum. 
This A'^min does not depend on the value of E°]!.lj^, and corresponds 
to a minimum fluence Emin. We have: 



En: 



= A^min 



e: 



peak 



Aff(-B°^Ik) 



(1) 



where AcS is the effective area of the detector. Inverting Eq. [T] 
we have a threshold curve _Epgak(^min): only for bursts with flu- 
ence larger than this curve we can derive E°JJak- For a constant 
AcH, bursts with larger Epcak have larger Fmin (since TVmin is the 
same), and this in principle might explain why only few faint bursts 
have large E°^^^. For a fixed Nmin and constant Acs, the limit- 
ing curve in the EpJjJk~Fltience plane is linear. This is indeed very 
similar to the shape of the limiting regions for BATSE, Swift, and 
BeppoSAX-GRBM reported in Fig. [5] 

In Fig.[5]we also show the limiting region for the spectral anal- 
ysis for SeppoSAX-GRBM. Although the response files of this in- 
strument has not been yet released, the GRBM and BATSE are sim- 
ilar scintillators and their effective area are similar in the ^20-800 
keV range. For this reason we can have an approximate idea of the 
limiting fluence by rescaling the limiting Epcak(-fmin) curves of 
BATSE for the ratio of the effective areas of BATSE and GRBM. 

We note that the "spectral threshold" is not affecting the sam- 
ple of BATSE and BeppoSAX bursts (solid open symbols in Fig. 
[Sj. It is therefore compelling to investigate if there exists bursts 
detected by BATSE or BeppoSAX which lie between the limiting 
threshold curves and the sample of GRBs with measured redshifts 
shown in Fig. [5] These bursts would not have a measured redshift 
and would not be affected by the selection effect (if any) interven- 
ing when asking that GRB has its redshift measured. In order to 
investigate these issues, however, it is necessary to have a com- 
plete GRB sample down to the limiting fluence represented by the 
curves of Fig. |5] and we leave this to a forthcoming paper (Nava et 
al. 2008). 

We can see that the "spectral threshold" is affecting the Swift 
bursts (filled red squares in Fig.|5}: due to the the limited BAT en- 
ergy range. Swift can only add bursts in the ~ 15-150 keV range of 
the _BpI^ak~Fltience plane. Note that there are 3 Swift bursts which 
are below the spectral thresholds defined above. One of these is 
GRB 060218 (open square in Fig.[5ll whose spectral peak energy 
has been found by combining the BAT and XRT Swift data (Cam- 
pana et al. 2006). The other two are GRB 050416A and GRB 
060526. For GRB 0504 16A we adopted the spectral results re- 
ported in Sato et al. (2007) which were derived through a fit with a 
Band model with a fixed low energy spectral index in order to ac- 
count for the softness of the burst (see also Sakamoto et al. 2006). 
The spectrum of this burst can be well fitted also with a single 
power law softer than F = 2. For GRB 060526 we adopted the 
spectral parameters reported in Schaefer et al. (2006). However, 



This is controlled by the pcorferf keyword of the batphasimerr tool 
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Figure 6. Simulation of the dependence of the slope of the B°''^. -Fluence 
correlation on the energy interval width \og(E^/E{) in S°Jj|j. used to 
compute it. Here Eu and Ei represent the boundary of the interval in 
^°cak' '^'^^ original correlation (with slope 0.39) and its 90% confidence 
level are shown by the solid and dotted orizontal lines, respectively. The 
solid line and the dashed lines represent the slope of the correlation obtained 
by selecting a sub-sample of bursts with peak energy comprised within a 
range of width log(i?u/£;i). 



our analysis of the BAT datqj were consistent also with a single 
power law with V ~ 1.7. 

Having shown that at least the low end of the correlation in 
the iJpJj^j.-Fluence is affected by a truncation effect, acting mainly 
on the Swift burst population, we proceed to analyze the correla- 
tion following the method proposed by LPMOO. This method, how- 
ever, can be applied only to the sub-samples of bursts for which the 
dominant truncation effect is known. We have shown that nor the 
trigger detection threshold neither the spectral analysis threshold 
of the different instruments that contributed to the heterogeneous 
GRB sample are affecting the GRBs with redshifts except for the 
Swift sample. We therefore can apply the correlation analysis on 
the Swift sub-sample taking into account its spectral analysis trun- 
cation effect. 

We calculate the Kendall's correlation coefficient with only 
those Swift pairs of bursts which lie above each other's spec- 
tral analysis thresholds. We obtained r = 0.1 (0.6(t significant). 
The significance has been computed with the boostrap method de- 
scribed in the Appendix of LPMOO. This result suggests that there 
is no iSpJ^li^-Fluence correlation in the Swift GRB sample consid- 
ered. However, we stress that the Swift bursts suffer from another 
strong truncation effect: the bursts peak energy can be computed 
only if it lies in the 15-150 keV energy range where the spectra are 
fitted (but see B07 for a different approach). This is clearly shown 
in Fig.|5]by the clustering of the Swift bursts (red squares) in this 
limited energy range. This energy range is smaller than the dis- 
persion of the i5°[^Jj.-Fluence correlation defined by the pre-Sw;// 
sample and this implies that it cannot be used to probe the exis- 
tence of the i?pgJi;-Fluence correlation. Indeed, the dispersion of 

* www.brera.inaf.it/utenti/gabriele/060526.html 



the Swift bursts along the fluence axis is much larger than the dis- 
persion along the E'^l^ axis only due to the limited energy range 
where the Swift burst peak energy is computed. To prove this we 
performed a simulation: assume that the correlation defined by the 
76 GRBs in the BpJjJi^ -Fluence plane is true and extract randomly 
a subsample of bursts with iJpeJit within the range [E\-E^]. We 
vary the width of this range and compute the correlation of the 
resulting subsample. The result is shown in Fig. [6] Note that the 
energy range in iSp^ak of Swift bursts (red squares in Fig. |5] ex- 
cluding 050416A, 060526 and 060218 - see above) corresponds to 
Eu/Ei ~ 5. With this small energy range we do not recover the 
true correlation. Only if Eu/Ei>25 we have a correlation consis- 
tent with the simulated one. Based on this result, we conclude that 
while the Swift bursts alone do not show the iSpj!^]^ -Fluence corre- 
lation shown by the pre-Swift bursts, this sample cannot be used to 
rule out the existence of the very same correlation due to its lim- 
ited range of Sp^ak- Considering instead all the other bursts we can 
conclude that the selection effects considered are not affecting their 
distribution in the i5pJ^Jj.-Fluence plane but still we cannot exlcude 
that other, unknown and still to be studied, selection effects could 
be present. In this paper we studied two of the most relevant selec- 
tion effects, i.e. the condition of detecting a burst and that (more 
constraining) of properly fitting its spectrum to derive E°^ly. and 
its fluence. 

Finally, there could be another selection effect responsible for 
the fact that the bursts of measured redshift have on average larger 
fluences. This selection effect could be associated to the require- 
ment of deriving their position with an accuracy high enough to 
start the ground based follow-up. Indeed, while in the Swift era the 
X-ray afterglow position is known with a typical few arcsec accu- 
racy, in the pre-5w(/f era the positioning of the X-ray transient was 
given by the WFC and WXM. In order to test this possibility we 
considered the sample of bursts detected by the WFC and GRBM 
on board BeppoSAX (Frontera 2004). In this case, as shown by 
the trigger threshold reported in Fig. [5] it was the GRBM trigger 
threshold to dominate over the WFC one. However, by compar- 
ing the distribution of the fluence of the sub-sample of BeppoSAX 
bursts with and without redshifts (Fig. |7), we see that there is no 
preference for bursts with redshifts to have larger fluences (the 
Kolmogorov-Smirnov probability test is 0.38). The bursts with red- 
shift are about half of those detected by both the GRBM and WFC 
onboard BeppoSAX. This is expected, due to the fraction of dark 
bursts, and the visibility constraints of the ground based optical 
telescopes. 

This result, regarding the BeppoSAX GRB sample, is puz- 
zling. In fact, in Fig.|5] the WFC limiting curves are a factor 100 on 
the left side of the distribution of the corresponding BeppoSAX 
bursts (also shown in Fig. |7]l. We would have expected that the 
GRBs detected by both the WFC and the GRBM, without requiring 
the redshift to be measured, had minimum fluences corresponding 
to the limit posed by the GRBM. Instead they all lie at larger flu- 
ences, and we did not find any reason to explain this puzzle. 

In Fig. [8] we show in the Speak — Siso plane, separately, the 
pie-Swift bursts and those whose spectral parameters have been de- 
termined only by Swift data. The best least square fit to the two sets 
of data yields a slope a — 0.51 ± 0.03 for the pve-Swift GRBs 
and a = 0.36 ± 0.05 for the Swift bursts. The two slopes become 
consistent if one performs a fit that weights for the errors on both 
quantities, as found by C07. However, for the Swift bursts of our 
sample, we have found the dominant selection effect. By account- 
ing for this truncation we find a very weak correlation in the rest 
frame (a = 0.1, Kendall's r = 0.12 at O.58o-). 
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Figure 7. Distribution of the observed fluence (2-700 keV) of bursts de- 
tected by both the GRBM and the WFC on board BeppoSAX (Frontera 
2004). The two histograms con'espond to GRBs with (hatched) and without 
(empty) redshifts. The two distribution Kolmogorov-Smimov test probabil- 
ity is 0.38. 
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Figure 8. Rest frame Bpcak — ^iso correlation. The Swift bursts whose 
peak energy has been derived from the analysis of the BAT data are shown 
with filled red squares (their fit - excluding GRB 060218 - is represented 
by the solid red line). The non-Swift bursts are shown by the open symbols 
and their fit by the long-dashed black line. 



5 SUMMARY AND CONCLUSIONS 

In this paper we have studied the Speak — Eiso correlation defined 
with the most updated sample (up to Sep. 2007) of GRBs with 
measured redshift and well defined spectral properties: 35 GRBs 
detected before the launch of the Swift satellite (Nov. 2004) and 41 
events detected since then. The latter sample is composed by 27 



events whose spectral properties were obtained through the analy- 
sis of the Swift-Bfiil spectra and 14 bursts (in most cases also de- 
tected by Swift) but with a spectrum detected by other satellites 
(Konus, Suzaku, RHESSI) with a larger spectral energy window. 
With this large GRB sample we have studied the possible evolution 
with redshift, the existence of other correlations among the spec- 
tral parameters or with redshift, the possible instrumental selection 
effects acting on this sample of 76 GRBs with measured redshift. 

Our main results show that in the rest frame Epak—Eiso plane 
the 76 GRBs define a strong correlation with no new outlier (with 
respect to the classical GRB 980425 and GRB 031203 (see e.g. 
Ghisellini et al. 2006). We find no correlation between the spectral 
parameters and the redshift nor between the peak energy and the 
spectral photon index. 

By dividing the GRB sample into redshift bins we compared 
the slope of the Speak correlation for each redshift bin. We do not 
find any evidence that this correlation evolves with redshift (as also 
found by B07) contrary to what has been found with a smaller GRB 
sample (Li 2007). 

The analysis of the observer frame Sp^ak-Fluence plane 
shows the existence of a strong correlation between the observed 
peak energy and the fluence. This correlation is of the form 
Speak oc(Fluence)''* ^rid GRBs of different redshift are evenly dis- 
tributed along it. The Sp^^i^ -Fluence plane is where the possible 
selection effects on the Speak correlation should be investigated. 
Two regions of this plane are defined by the Speak~Fluence corre- 
lation: bursts with large fluence and small/intermediate peak energy 
should not exist (or be very rare) as nothing prevents their detection, 
while bursts of intermediate/large Speak ^rid small fluence could 
be affected by some instrumental selection effect. We consider two 
of these: (a) the minimum flux to trigger a burst and (b) the min- 
imum fluence required to fit its spectrum and constrain its peak 
energy. We have used the threshold curves derived by Band (2003) 
to represent the first selection effect in the Speak-Fluence plane. 
We have converted these curves into the fluence-Speak plane, as- 
suming a range of burst durations. We found that the trigger thresh- 
old does not affect the distribution of bursts with known redshifts 
in the SpeJ^-Fluence plane. The second selection effect that we 
considered (to our knowledge this is the first time that such a se- 
lection effect is considered for the GRB detectors) is related to the 
requirement to have enough photons in the spectrum to fit it with 
a curved spectral model (either the Band model - Band et al. 1996 
- or a cutoff-power law model) and constrain the peak energy of 
the burst. To this aim we performed detailed spectral simulations 
which account for the detector response and the typical GRB back- 
ground. We also accounted for some parameters (burst duration and 
GRB off-axis position) which contribute to determine these thresh- 
old curves. Our results (shown in Fig.|5l( demonstrate that the latter 
selection effect dominates the first in the SpeJ]j-Fluence plane as 
shown in Fig. [5] Our results show that the pie-Swift GRB sample, 
containing a fraction of bursts detected by BeppoSAX and BATSE 
is not affected by the corresponding limiting curves, while the Swift 
sample is. Indeed, in the case of Swift most of the GRB whose spec- 
tra are well fitted by a cutoff-power law model have an S°e|]j which 
is in the range 15-150 keV. The correlation in the Speak-Fluence 
plane defined by the Swift sample is flat and very weak compared 
to that defined with all the other bursts. However, we cautiously 
note that the Swift sample is distributed in a very narrow range in 
Speak, smaller than the scatter of the Speak-Fluence correlation 
as defined by all the pie-Swift bursts. As a consequence, although 
Swift bursts do not show a Speak-Fluence correlation, we cannot 
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rule out the existence of the very same correlation defined by the 
non-Swift bursts, extending over 2 orders of magnitudes in both 
£;°^:kandFluence. 

Our results do not exclude that other selection effects affect the 
^peak~Fluence plane, and in particular the non-Swift bursts. The 
lack of many bursts with known redshift with intermediate/large 
Spcak^nd small fluence (i.e. between the present sample of GRBs 
and the BATSE limiting curves of Fig. [5j could be due to the ad- 
ditional presence of a still not understood selection effect, different 
from those considered in the present paper (see also Nakar & Pi- 
ran 2005). The first step to investigate this possibility is to search if 
there exist GRBs which populate the region on the left side of the 
£°j!|i5 -Fluence correlation. This issue cannot be solved by includ- 
ing, in the i^pj^lj^ -Fluence plot, the BATSE bursts already spectrally 
analyzed by Kaneko et al. (2006), since they all have large fluences. 
Furthermore, the sample of Yonetoku et al. (2004) could be biased 
by their requirement of having a pseudo-redshift not exceeding 12, 
as measured by the _Bpeak-luminosity relation. Therefore we need a 
new, representative sample of BATSE bursts with measured Sp^Jk 
and fluences close to the limiting fluence curve (see again Fig.[5]l. 
This is what we will present in a forthcoming paper (Nava et al. 
2008, in preparation). We anticipate here that these bursts exist and 
therefore the conclusion that can be drawn from the present work 
is that other selection effects, different from the the ORB trigger 
threshold or the spectral threshold considered in the present paper, 
are very likely affecting the iJpJjJi^-Fluence correlation and, as a 
consequence, the rest frame E-pcak ~ Eiso correlation. 
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